Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

Molecular flow and wall collision age distributions

The use of storage cells has become a standard technique for internal gas targets in conjunction with high energy storage rings. In case of spin-polarized hydrogen and deuterium gas targets the interaction of the injected atoms with the walls of the storage cell can lead to depolarization and recombination. Thus the number of wall collisions of the atoms in the target gas is important for modeling the processes of spin relaxation and recombination. It is shown in this article that the diffusion process of rarefied gases in long tubes or storage cells can be described with the help of the one-dimensional diffusion equation. Mathematical methods are presented that allow one to calculate collision age distributions (CAD) and their moments analytically. These methods provide a better understanding of the different aspects of diffusion than Monte Carlo calculations. Additionally it is shown that measurements of the atomic density or polarization of a gas sample taken from the center of the tube allow one to determine the possible range of the corresponding density weighted average values along the tube. The calculations are applied to the storage cell geometry of the HERMES internal polarized hydrogen and deuterium gas target. Received 9 July 2001 and Received in final form 18 September 2001